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ABSTRACT 

Aims. We study the connection between the chromospheric and photospheric behaviour of the active late-type star FK Comae. 
Methods. We use spot temperature modelling, light curve inversion based on narrow- and wide-band photometric measurements, 
Ho- observations from 1997-2010, and Doppler maps from 2004-2010 to compare the behaviour of chromospheric and photospheric 
features. 

Results. Investigating low-resolution Hu spectra we find that the changes in the chromosphere seem to happen mainly on a time 
scale longer than a few hours, but shorter variations were also observed. According to the Ho measurements prominences are often 
found in the chromosphere that reach to more than a stellar radius and are stable for weeks, and which seem to be often, but not 
every time connected with dark photospheric spots. The rotational modulation of the Ho- emission seems to typically be anticorrelated 
with the light curve, but we did not find convincing evidence of a clear connection in the long-term trends of the Hu emission and 
the brightness of the star. In addition, FKCom seems to be in an unusually quiet state in 2009-2010 with very little chromospheric 
activity and low spot contrast, that might indicate the long-term decrease of activity. 

Key words. stars: activity - stars: magnetic field - stars: atmospheres - stars: chromospheres - stars: late-type - 
stars:individual:FK Com 


1. Introduction 


FK Comae is the eponymous member of a group of active stars 
( jBopp & Rucinski||1981| ), where the members are fast-rotating 
G-K giants showing chromospheric and coronal activity signs 
similar to the RS CVn-type binaries. Their typical v sin i is of or¬ 
der of 100km s^ 1 , but no periodic radial velocity variations have 
been detected, indicating that these stars are single, or have a 


companion of very low mass. Studying the Na D line, McCarthy 
& Ramsey! (1984) gave 5 km s 1 as an upper limit of the radial 
velocity variations caused by a possible companion, correspond¬ 
ing to a mass of 0.054 Ms im . Huenemoerder et al. ( 1993j ) gave 


an upper limit of the semi-amplitude of 3 kms -1 , which makes 
the presence of the companion even more unlikely. Photomet 
ric observations showed periodic, quasi-sinusoidal changes of 
AV = O'f 1 - 0™3 interpreted as starspots by Bopp & Rucinski 


( 1981[ >. Bopp & Rucinski ( |1981| also suggested, based on the 
rapid rotation and the lack of radial velocity variations, that these 
stars are coalesced W UMa-type systems following the scenario 
proposed by Webbink (1976). 


* Based on the observations obtained at the Nordic Optical Tele¬ 
scope, Observatorio Roque de los Muchachos, La Palma, Canary Is¬ 
lands, Spain; Kitt Peak National Observatory, USA 


FK Com is the most active, and most thoroughly studied star 
of this group. The star shows signs of activity from photosphere 
to corona (|Merrill|1948||Chugainov||1966||Ramsey et al.j | 1981 


Dorren et al. 1984; Oliveira & Foing| [1999[ |Ayres et al. 2006 

Drake et al. 2008). FK Com was also the first object where the 


so-called flip-flop phenomenon was observed ( |Jetsu et aT7| 1993 
|1994| ): the star shows active longitudes separated by -180°, and 
the dominance between these regions is exchanged every few 


years. Olah et al. (2006) found that both flip-flops and phase 


jumps can be found on FK Com. They concluded, that these are 
two different, although possibly connected phenomena. |Hack-| 
man et al. (|2013|> also showed that flip-flops are not a single 


phenomenon: they can occur gradually (through differential ro¬ 
tation) or suddenly, the phase shift varies and also as does the 
stability of the new primary spot region. 


Ramsey et al. ( 198 lj) described the broad, double-peaked 
Her profiles of FK Com with an “excretion” disk around the 
star, which itself has a small contribution to the Her emission. 
Huenemoerder et al. (1993}) presented photometric and spectro¬ 


scopic data and concluded that FK Com is a single star with ex 
tended matter around it. From long-term Her observations jWelty 


et al. (jl993)> found evidence of co-rotating emitting material 


around FK Com within two stellar radii with lifetimes longer 
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Table 1. List of Ha observations 


High-resolution spectra 


Date 

Telescope 

# of spectra 


1997/04/07-04/16 

KPNO 

10 


1998/03/10-03/17 

NOT 

7 


1998/04/01-04/22 

KPNO 

19 


1998/07/02-07/15 

NOT 

38 


1999/03/03-03/05 

NOT 

6 


1999/05/24-06/04 

NOT 

19 


1999/07/23-08/03 

NOT 

12 


2000/03/29-04/08 

KPNO 

10 


2000/04/26-05/05 

KPNO 

10 


2000/08/07-08/18 

NOT 

12 


2001/05/06-05/09 

NOT 

6 


2001/06/03-06/12 

NOT 

17 


2002/08/20-08/29 

NOT 

7 


2003/06/03-06/22 

NOT 

13 


2004/01/30-02/11 

NOT 

10 


2004/07/25-08/08 

NOT 

10 


2005/07/15-07/23 

NOT 

7 


2007/07/18-07/29 

NOT 

12 


2009/08/26-09/06 

NOT 

11 


2009/12/27-01/04 

NOT 

18 


Low-resolution spectra 


Date 

Telescope 

# of spectra 

Phase 

1999/02/06 

NOT 

11 

0.35-0.40 

1999/02/07 

NOT 

43 

0.75-0.78 


than a few weeks. Oliveira & Foing (1999) modelled highly vari¬ 
able Balmer line profiles and confirmed the complex circumstel- 
lar structures. [Kjurkchieva & Marchevj ( j2005j ) presented high- 
resolution Ha observations spanning over two years. The authors 
suggested a model where FK Com is a binary system with a low- 
mass secondary star illuminating a circumstellar accretion disk. 
The model contains additional sources of absorption and emis¬ 
sion at phases 0.0 and 0.5, respectively. 

In this paper, the seventh in a series on FK Comae, we anal¬ 
yse Ha data spanning over a decade from 1997 to 2007. The 
chromospheric activity patterns from the Ha data are compared, 
in most epochs, to the contemporaneous photospheric spot maps 
obtained using the Doppler imaging technique. 


2. Observations 


Most of the spectroscopic observations presented here were ob¬ 
tained at the 2.56 m Nordic Optical Telescope (NOT) during 13 
observing runs using the SOFIN high resolution spectrograph 
and one run of low resolution spectroscopy using ALFOSC. In 
addition, three datasets of high resolution spectra were obtained 
at the Kitt Peak National Observatory (KPNO) in Arizona. All 
the runs from KPNO and 10 runs with the SOFIN at NOT (1998- 
2003 and 2004 February) were already used in Olah et al7| ( 2006[ >, 
|Ayres et ah] ( |2006| ) and |Korhonen et al.| ( |2007 ! for investigating 
the photospheric activity. Thus, here only the unpublished ob¬ 
servations from SOFIN and ALFOSC are described in detail. 
In addition to these spectroscopic observations, also previously 
published photometric observations from Korhonen et al. ([2001) 
and |Hackman et alT| ( |2013[ ) are used. 

The phases for all the observations presented here were cal¬ 
culated using the ephemeris obtained from 25 years of photo¬ 
metric observations, HJD = 2,439,252.895 + (2‘?4002466 ± 
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0^0000056 )E, referring to a photometric minimum calculated by 
IJetsuetal. 1 ( 1993 , 1994). 

2.1. High resolution spectroscopy 

High resolution NOT observations for five epochs (2004 July, 
2005 July, 2007 July, 2009 August and 2010 January) were car¬ 
ried out using SOFIN echelle spectrograph with the low reso¬ 
lution camera, giving resolving power (A/AA) of 27,000. The 
echelle spectra consisted of typically 45 orders and were cen¬ 
tred at 5200 A. The signal-to-noise ratio in the Ha order varied 
between 90 and 310, typical value being around 200. The data 
were reduced using the 4A reduction package ( |Ilyin|2000| . Ta¬ 
ble [T| gives the summary of the high resolution spectra used in 
this paper, including the previously published datasets. 

2.2. Low resolution spectroscopy 

For investigating the possible rapid variations in the Ha profiles, 
low resolution spectra with high cadence were obtained using 
ALFOSC at the NOT during 1999 February 26-27. We used the 
0.4 arcsec slit and the echelle grism #9 together with the Ha filter 
#22. With this configuration we selected only the order contain¬ 
ing Ha, without using a cross disperser. The resolving power 
obtained with this configurations was approximately 4,000. 

The observations were done in sets of two observations with 
60 sec exposure times, and taking Halogen flat and Helium- 
Neon arc observations before and after each pair. Due to se¬ 
lecting only a narrow spectral region only a couple of arc lines 
were seen. Thus, no accurate wavelength calibration could be 
carried out. An approximate wavelength solution was obtained 
from the theoretical calculations for the grism #9, and the solu¬ 
tion was fine-tuned for each observation by cross-correlating the 
arc spectra and introducing the shifts obtained to corresponding 
observation. As we were only interested in the rapid variability 
of the profiles this procedure gave a wavelength solution that was 
sufficient for our purpose. All the observations were reduced us¬ 
ing the Image Reduction and Analysis Facility (IRAF) which is 
distributed by KPNO/NOAO. 

2.3. Photometry 

In this paper beside the AV observation presented in 
( |2QQ9[ > and |Hackman et aI7| ( |20131 >, we use the all the 
in different other colours (U, B, Ic, b, y) parallel with the V mea¬ 
surements of |Korhonen et al.| ( [2001[|2007j l, to estimate the tem¬ 
peratures of active regions in FK Com. 

3. Results 

3.1. Long term chromospheric activity 

A sample of the Ha spectra is plotted on Fig.[T| For each observ¬ 
ing run 5-10 spectra that show well the behaviour of the region 
were chosen and plotted. As can be seen, the profiles are highly 
variable, but at each epoch display a broad double-peaked emis¬ 
sion which is indicative of a circumstellar disc. 

The complex Ha profile shape of FK Com indicates many 
different components: circumstellar disc, plages and promi¬ 
nences. To study just the contribution from the stellar component 
(plages and prominences) the disc contribution is removed. To 
facilitate this, we subtracted the spectrum that showed the least 
amount of stellar Ha emission and seemed to only show the cir¬ 
cumstellar disc component (a NOT spectrum from 2000/08/10) 


Olah et al. 
data taken 
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Fig. 1. Samples of high-resolution Her profiles from different observing runs. 


having a symmetric Ho- profile with no other emission feature 
than the two circumstellar disc-peaks. However, it is possible 
that this lowest intensity spectrum is just a local activity min¬ 
imum in the decade-long observations and even this one con¬ 
tains some low-level plage- or prominence-like chromospheric 
activity. Using these residual spectra, we have plotted the dy¬ 
namic spectra on Fig. [2] for each observing run. The plots show 
the excess emission at different velocities (x-axis) and rotational 
phased (y-axis, times of observations marked with plus-signs). 
The dashed white line denotes the location of the stellar disc 
centre and the dotted lines the edges of the stellar disc. Colour 
coding (black-red-green-blue) shows the intensity of the Ho- re¬ 
gion. Crosses mark the phases of the observations, for phases 
where there is no data the plot is interpolated to the closest 
measurements. More information on how these dynamical spec¬ 
tra have been calculated can be found from Korh onen et al.l 
(2009) . On the whole, most of the excess emission is seen out¬ 
side the stellar disc, implying material located higher in the stel¬ 
lar atmosphere. Sometimes, excess emission is also seen on the 
stellar disc, implying plage-like region in the chromosphere of 
FK Com. The wave-like features seen on these plots are proba¬ 
bly the result of prominences rotating well outside the stellar ra¬ 
dius. This feature can be best seen at 1998.27, 1998.30, 1999.41, 
2000.33, 2001.35, 2004.10, 2004.58, and 2007.56. The shape of 
the waves is continuously changing, indicating that these promi¬ 
nences are constantly evolving. It is possible, that some of the 
prominence-like features is related to the circumstellar disc: the 
disc could get heated unevenly, for example by stellar flares. 
This could cause some transient features in the disc emission. 
We cannot completely rule out that some of the emission seen 


would come from the disc, but the activity of the star itself is 
more likely the origin (the star is very active, as shown by all 
indicators). 


To check the long-term behaviour of the Her line, and give a 
more quantitative description, we fitted Gaussians to the residual 
spectra for measuring the changes of the whole Her region, and a 
simple model consisting of two Gaussians to the original spectra 
to describe the violet and red peaks of the Her region individ¬ 
ually. In Fig. [3] we plotted the AV light curve from Olah et al. 


(2009) and Hackman et al.| ( |201 3| > together with the strength of 
the Her emission measured as the area of the fitted Gaussians 
(both from the residuals and the individually fitted violet and 
red peaks). The strength of the line decreases until -2002-2003, 
then increases in all cases, indicating a change of chromospheric 
activity on a time scale of a few years. 


3.2. Fast variability in chromospheric activity 

Along with the long-term variations, we have also searched for 
minute scale variations using lower resolution spectra obtained 
in 1999 February (see Table [TJ. The spectra from the two nights 
are plotted in Fig. [4] The Her profile of FKCom does not show 
relevant changes during the time span of about 2-3 hours of each 
night. The observations from the first night show more small 
scale variations. To check if these changes are real or caused 
by the lower signal-to-noise ratio of these observations (average 
S/N was 88 on the first night and 118 on the second), we av¬ 
eraged the spectra in 30-minute windows, and plotted the result 
in Fig. [4] These averages indicate small variations in the Her re¬ 
gion on the time scale of hours, especially during the first night 
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Fig. 2. Dynamic Htt residual spectra of FK Comae. Crosses on the left show the phases of individual spectra. Thick dashed lines mark the 0 km s 
velocity shift (i.e., the centre of the stellar disk), the thin lines correspond to the edges of the disk. 
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Fig. 3. From top to bottom: AV photometry of FK Comae from |Qlah 
|et al.| < |2009| > and |Hackman et al.| l |2013} (the latter after 2008, with 
darker colour), Hu line areas from Gaussian fits of the residual spec¬ 
tra, and the area change of the Ho: violet (V) and red (R) peaks (V peak 
shifted arbitrarily by 20 A for clarity). 



Wavelength (A) 


Fig. 4. Low-resolution Ho from 1999 February. Red continuous line 
shows the first, blue dashed line shows the second night. The first night 
covers ~ 3, the second night ~ 2 hours. The same spectra, averaged in 
30 min windows, are plotted above the original, with colours brighten¬ 
ing with time. 


which is exhibiting higher activity levels. Anyhow, there is more 
prominent difference between the two nights. Thus we can con¬ 
clude that, at least during the two nights of our observing run, 
the chromospheric structure of the star shows significant night- 
to-night variability, and much weaker variations on smaller time 
scales. Note also that the low resolution of the observations can 
mask some smaller scale variability. 


3.3. Photospheric spots from Doppler imaging 


For the five previously unpublished high resolution spectral 
datasets we also obtain surface temperature maps using Doppler 
imaging techniques. For this the Tikhonov Regularisation in¬ 
version code INVERS7PD, written by N. Piskunov (Piskunov 
et al.| 1990 > and modified by T. Hackman (Hackman et al. |2001| 
was used. The adopted stellar parameters and additional details 
on the inversion technique and line-profile calculations can be 
found, e.g., in|Korhonen et al.|(|2000 12007)1. The obtained tem- 


perature maps are shown in Fig.pland discussed in the following. 

In 2004 July observations from 10 rotational phases of FK 
Com were obtained. The phase coverage of these observations 
is good, and the largest phase gap is 0.17, spanning the phases 
0.36-0.52. The map obtained from the data shows two high lat¬ 
itude, 47-79°, spot concentrations around the phases 0.25 and 
0.85. A lower latitude active region is seen around the phase 
0.15 spanning the latitudes 20-29°. The average spot latitude ob¬ 
tained from this map is 54.9 + 3.4° according to the output of the 
inversion code. The average spot latitudes were determined from 
the Doppler images using areas that were cooler than the limiting 
temperature of T\{ m = 4550 K (77 n = 5000 K unspotted surface 
temperature assumed), its error was estimated by the standard 
error of the weighted mean. The area at individual latitude strips 
were measured, and in the averaging the spot areas at each lati¬ 
tude were used as weights. The two spots located at the phases 
0.1-0.3 are approximately 1200 K cooler than the unspotted sur¬ 
face, and the one at the phase 0.8 is 800 K cooler. The coolest 
temperatures in the spot group around the phases 0.1-0.3 are up 
to 1500 K cooler than that of the unspotted surface. 

The SOFIN observations obtained in 2005 July consist of 
eight different phases, with the largest gap spanning the phases 
0.17-0.45. Not many cool regions are seen in this epoch. Only 
one spot, which is approximately 500 K cooler than the unspot¬ 
ted surface, is present at phases 0.7-0.8 spanning the latitude 
range 56-70°. The average latitude of this region is 61.1 + 0.7°. 
The coolest temperature in this spot region is about 600 K cooler 
than that of the unspotted surface. 

The phase coverage of the 12 spectra obtained in 2007 July is 
very good and there are no phase gaps larger than 0.1. Also this 
map shows, similarly to the 2004 July one, two spot concentra¬ 
tions on the surface. The main spot group is approximately 600 K 
cooler than the unspotted surface, with the coolest temperatures 
going down to 800 K cooler than the unspotted surface. This spot 
is located at the phases 0.7-0.8 and spans the latitudes 56-74°. 
The second spot group has similar temperatures and it is located 
at phases 0.1-0.2 and spans latitudes 61-70°. The average spot 
latitude is 65.3 + 0.8°. 

The observations carried out in 2009 August consist of 11 
spectra, which are relatively evenly distributed over the stellar 
rotation. The largest phase gap of 0.17 in phase is located at the 
phases 0.51-0.68. Similarly to the July 2005 map there are no 
prominent cool spots at this epoch. The coolest regions, which 
occur around phases 0.25-0.5, are only approximately 500 K 
less than the temperature of the unspotted surface. The latitude 
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of the spots is also relatively low compared to what is typically 
seen in FK Com, 7-34° with and average latitude of 18.5 ± 1.5°. 

In 2010 January 18 spectra were obtained on three nights 
during a four night time period. Therefore, even though there are 
many phases, they are grouped together and there are large phase 
gaps in the dataset. The largest gap occurs between the phases 
0.59 and 0.93. As in the August 2009 and July 2005 map the 
cool spots are of very low contrast and the difference between 
them and the unspotted surface is only approximately 550 K. 
Two main spot groups can be seen on the surface, one spanning 
phases 0.15-0.25 and another one around phase 0.8. As is the 
case for August 2009, also in January 2010 the spots are located 
at low latitudes, 2-25° with and average latitude of 12.4 + 1.2°. 


4. Discussion 


4. 1. Photospheric activity 

The evolution of the photospheric activity of FK Com in 2004— 
2010 can be studied based on the five new Doppler images pre¬ 
sented in this paper. 

The first map, July 2004, has the highest spot contrast of 
all the new maps presented here. Second highest spot contrast 
is seen in the July 2007 map, where a spot contrast typical of 
FK Com of approximately 800 K is seen. July 2005, August 
2009, and January 2010 maps, on the other hand, have similarly 
low spot contrasts with the coolest spots being only 500-600 K 
cooler than the unspotted surface. Often this kind of low con¬ 
trast maps are produced if the data quality is bad, but here espe¬ 
cially the July 2007 and August 2009 maps are of good quality, 
and particularly have very good phase coverages. Therefore, it 
is clear that the temperature of the spots has drastically changed 
between July 2004 and July 2005, and that also in the latter years 
FK Com was still in the low activity state. Note, that hot features 
in the maps at the same longitude as the prominent cool spots are 
most likely artefacts caused by changing continuum levels with 
changing temperature. If the hot features are well separate from 
the cool spots they could be real (see Section 4.4 1 . 


All the maps show the spots at similar phases, around 0.2 
and 0.8. All the spots are of very low contrast in the July 2005, 
August 2009 and January 2010 maps and it is debatable if they 
can even be considered active regions. The spot group around 
phase 0.8-0.9 is seen in all the maps, except August 2009, but 
in the July 2004 map it is the secondary spot group and the 
main group in the July 2005 and the July 2007 maps. This could 
indicate a possible flip-flop event between July 2004 and July 
2005, an event which is supported by the July 2007 map. On the 
other hand, the two later maps, August 2009 and January 2010, 
again show 0.2 as the more active phase. This would imply yet 
another flip-flop, which would have taken place between 2007 
and 2009. Korhonen et ak ( |2002[ ) reported that flip-flop events 
occur on average every 3.2 years, which would mean that two 
flip-flops could in principle have happened in the 5.5 years the 
observations cover (although the flip-flops are not regular phe¬ 
nomena, see e,g. |Qlah et al.|2006[|Hackman et al.|2013| l. Still, as 
is also concluded by Korhonen et al. 1 2007) 1, flip-flops are diffi¬ 
cult to detect from sparse Doppler imaging alone, and are better 
seen based on photometric observations with denser coverage. 
For this reason we also carried out light curve inversion, where 
inversion techniques are applied to broad-band photometry and 
maps of the fraction of each pixel on stellar surface covered by 
spots, so-called spot filling factor, are calculated. More on the 
technique can be found, e.g., from|Lanza et al.|(jl998)>. 
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The light curve inversion results are shown in Fig. [6] The 
data were divided into 21 subsets based on the stability of the 
light curve. Continuous phase shifts of the spotted longitudes can 
be followed in 6 subsequent maps between the JD 2453843.94- 
2454254.98 maps with one seasonal gap of 150 days in the data. 
Another continuous shifting of the spots is observed on 4 consec¬ 
utive maps between 2454520.78 and 2454623.75, corresponding 
to data without long interruptions. Drastic difference is seen be¬ 
tween the maps of 2454254.98 and 2454464.92, and 2454623.75 
and 2454835.65, with seasonal gaps in the data of about 150 
days. These cases, in lack of more continuous sequences, the 
changes can also be explained by continuous phase shifts, but 
could not be followed. However, a substantial change of the ac¬ 
tivity pattern, although shown by low amplitude light curves, 
happened between 2454867.69 and 2454934.47 maps, with a 
small gap of 28 days in the data, i.e., sometime in 2009 March, 
which can be considered as a flip-flop-like event (see the Discus¬ 
sion later). 

The spot filling factor maps can be compared with the 
Doppler images. The spot filling factor map closest to the July 
2005 is the one from Julian Date 2453473.08 (second map in 
the plot). The location of the main spot measured from the spot 
filling factor map is 0.73 in phase, similarly to what is seen in 
the Doppler image. The July 2007 Doppler image is close to the 
spot filling factor map from Julian Date 2454254.98. Again the 
location of the main spot in both maps is very similar. For Au¬ 
gust 2009 there is no contemporaneous spot filling factor map, 
but the map of 2454970.75 is obtained few months before it. 
The Doppler image of August 2009 shows very little spots, but 
they are located at same phases as the main activity seen in the 
spot filling factor map. Naturally, the activity could have gotten 
weaker within the approximately two months separating the spot 
filling factor map and the Doppler image. The spot filling factor 
map of Julian Date 2455247.87 is contemporaneous to the Jan¬ 
uary 2010 Doppler image. The spot filling factor map, like the 
Doppler image, shows main activity around the phase 0.25, and 
the surface could be covered by relatively low contrast spots. 

On the whole the spot filling factor maps obtained from the 
photometry correspond very well to the Doppler images that are 
based on high resolution spectra. The possible flip-flop events 
seen in the Doppler images are also seen in the spot filling factor 
maps, but suspiciously most of them occur during the long gap 
when the star is not observable. 

4.2. Correlating chromospheric and photospheric activity 
regions 

The low-resolution Her spectra from 1999 were observed be¬ 
tween phases 0.35-0.40 on the first night and 0.75-0.78 on the 
second night. Although no changes could be observed during 
the nights, there is considerable difference between the spectra 
of the two nights. On the first night the Her emission is stronger 
than on the second one. |Korhonen et al,| ( |2002| > published light 
curves and light curve inversions of FK Com from 1979 to 2001. 
The light curve from 1999.21 is obtained just a few days after 
the low-resolution Her spectra. In the light curve two spots can 
be seen: one at phases between 0.3-0.6, and one between 0.9- 
0.1. The Doppler map from 1999 April in Korhonen et al. ( |2007| > 
also shows that the most spotted state is around phase 0.3, where 
the star showed high-latitude spots. Thus the spectra of the first 
night show a spotted phase, while those from the second night 
show a less-spotted phase of FK Com. This can indicate a con¬ 
nection between the photospheric and chromospheric activity, 
as observed on the Sun and also on other stars: e.g. on the BY 
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Fig. 5. The surface temperature maps for 2004 July, 2005 July, 2007 July, 2009 August and 2010 January datasets. The lines above the maps mark 
the phases of the observations. 


Dra-type EY Dra (e.g., Korhonen et al. 2010 1 on the T Tauri-type 
TWA 6 qSkelly et al.|2008| >; on the K-dwarf LQ Hya ( |Frasca et al. 


2008); and on RS CVn-type binaries, e.g. RT Lacertae (Frasca 


et al.|2002] ); and on the W UMa-type VW Cep ( jFrasca et al. 


19961. The Her emission is however not necessarily connected 


to the photospheric spots: on SAO 51891 no modulation of Her 
was observed, although changes caused by chromospheric inho¬ 
mogeneities were clearly detected in the V light curve, and in the 
Ca ii HK, IRT and He lines. In this case the Her modulation might 
be possibly hidden by other activity signatures, e.g. microflares 
(see |Biazzo et al.|2009| l. 

It is interesting to note, that the spectra taken at the spotted 
phase of the star have asymmetric shapes, with a stronger peak in 
the red side of the Her line. |Kjurkchieva & Marc hey] ( |2005[ i found 
that the violet peak of the Her spectra shows a quasi-sinusoidal 
change with a maximum at phase 0.75 and a minimum at 0.25. 
The red peak was in maximum at phase 0.25 and in minimum at 
0.75, but the changes were not so significant as at the violet one. 
During this time, the main spotted area could be seen between 
phases 0.5-1.0 (see |Korhonen et al.|2004j) , i.e. at the light curve 
minimum the violet peak was in maximum. 


black arrows outside the circle - by visual inspection of the dy¬ 
namic Ha spectra shown in Fig. [2] selecting the phase when the 
Doppler shift of the prominence is zero. Sudden increases (prob¬ 
ably caused by flares) of Ha emission are also shown, marked by 
red arcs or dots. 

In most cases, when the wave-like behaviour of the Ha 
region was observed, the prominence appeared at a spotted 
phase, except at 2004.10. In many cases, the phases of the 
prominences coincide with the increased level of Ha activity; 
in 1998.30, 2000.33, 2004.10, 2004.58. At 1998.30, 1999.58, 
2001.35, 2004.58, and 2007.56 the prominence is very close 
in phase to a dark spotted feature seen in the Doppler-map. 
These can be an indication that there is a connection between 
the prominence and the photometric spots as found by Huen- 
|emoerder et al.| ( p~993| >, who also observed in 1989 that the dark 
photometric spots and Ha excess coincide in phase. 

A detailed description of the possible connection between 
the different activity features is given in Appendix [A] 


4.4. General level of activity 


4.3. Prominences 


In some cases in Fig.[2]noticeable wave-like features can be seen 
in the residual spectra. These are probably caused by promi¬ 
nences in the outer stellar atmosphere. According to the maxi¬ 
mum amplitude of these waves the emission mostly originates 
from about one stellar radius above the surface. In some cases 
of Fig. [2] (1999.58, and possibly in 1999.41, 2001.43) multiple 
waves can be seen, possibly caused by prominence arcs of differ¬ 
ent size and location. We compared these spectra to the photo- 


spheric features seen in the nearby light curve inversions (Korho- 
nen et al.f2 002 and this paper) and Doppler-images (Korhonen 
et al.|2007 and this paper in Fig. [5]). In Fig. [7] the phases of the 
spots, flares, and prominences can be seen. We marked photo- 
spheric spot locations both from light curve inversions and from 
Doppler images (with black arcs and green arcs/dots). We deter¬ 
mined the approximate phases of the prominences - marked by 


To measure the strength of the Ha line, we fitted Gaussians to the 
residual of the Ha region after subtracting the spectrum showing 
the lowest activity, thus removing the effects of a possible disk 
present around the star. 

Anti-correlated behaviour was observed between the V light 
curve and Ha by Huenemoerder et al. ( jl993) >, and this was the 
case in our low-resolution spectroscopic observations, too. If we 
compare phased light curves and the area of the residual Gaus¬ 
sians (see Fig. [8J, we can find similar trends, although the rela¬ 
tion here is much less clear, and was found only in some cases, 
while positive correlations and no Ha modulation have also been 
observed. One reason for this could be that the Ha emission 
is originating from both plage regions and prominences, and 
these features can be observed better at different phases. Namely, 
plage-like structures are better seen on the stellar disc, while 
prominences show themselves better when seen off-limb (see 
e.g. |Hall & Ramcey|[l992)l. Note, that strict anticorrelation be- 
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Fig. 6. The spot filling factor maps for 2005 January - 2010 April. In light curve plots the crosses denote the observations and the solid line the 
inversion result. In the spot filling factor maps the darker colour means larger spot filling-factor. In both plots the x-axis is longitude in degrees and 
the y-axis in the light curve plots is V magnitude and in the spot filling factor maps latitude in degrees. The mean Julian Date (HJD-2400000) of 
the observations for each set is given on top of the light curve plots. 


tween the V light curve and Her is observed only when the spots 
and the H<r emitting regions are co-spatial. 

Comparing the H<r intensities and the AC light curve (see 
Fig. [3J we find, that they have a very similar trend between 2000 
and 2006, although this connection cannot be seen in the ear¬ 
lier and later observations. The long-term positive correlation of 
the light curve and the Her line (i.e., with a stronger chromo¬ 
spheric activity as the star gets brighter) could be explained by 
facula-dominated activity, i.e., that most of the optical brightness 
changes are caused by hot spots - facular regions, at least on the 


long-term base. The magnitude-colour diagrams of Korhonen 
et al.] ( |2001] > showed that activity of FK Com is dominated by 
spottedness and not photospheric facular regions: the star gets 
redder as it becomes fainter (see also |Messina[2008T ). This find¬ 
ing can be confirmed by reconstructing these magnitude-colour 
diagrams using also the more recent photometric data of [Hack-] 
man et al. ( 2013| ) see Fig.[9]and Fig. [TO] These diagrams though, 
show the cumulative effect of the rotational modulation and the 
long-term changes which can hardly be separated. 


In Fig. 11a two-period fit to the V data is drawn in the upper 


panel, with lengths of 6.8 + 0.5 and 4.1 +0.1 years (the errors of 
the periods were estimated by the frequencies of the Fourier- 


spectrum 0 '."0005 below its maximum which corresponds to 
10% of the precision of the data used, cf. Press et afj 1992 ). 
A fit with the same periods is plotted on Ha area dataset, and de¬ 
scribes it reasonably well, apart from the first data (1997) with 
some phase shift. These two periods (rather time-scales) were 
found from a simple Fourier analysis (see online Fig. [T2) , and 
serve only to show the long-term changes in Her relative to V. 
However, these time-scales are verified with a time-frequency 
analysis of the 32 years long V light curve taken between 1979- 
2011, with the STFT (Short-time Fourier transform). The results 
are shown in the o nline Fig.[l3l For the details of this method and 
its application see |Kollath & O lah (2009l and |Qlah et al.| p009). 
After JD 2451000, i.e., around 1999 the long-term, mean bright¬ 
ness of the star shows an abrupt decrease. The cycle period is in¬ 
creasing from about 5.5 years in 1982 (2445000) to the present 
8.3 years. In 1999 the cycle period was about 6.6 years and a 
constant other cycle showed up with an approximate length of 
4.2 years, lasting till the present. This other cycle showed up af¬ 
ter the flip-flop event in 1997-1998. Whether any change in the 
cycle periods appear after the recent (2009) event remains to be 
seen in 15-20 years. 


Article number, page 8 of 14 





































































K. Vida et al.: Study of FK Comae Berenices. VII. 





Fig. 7. Possible connection between the dark photospheric spots and the 
chromospheric features. Approximate location of starspots from light 
curve inversion jKorhonen et al.|[2002fr and from this paper (see Fig. 
6|> are marked w ith black arcs. Stronger spots from Doppler-images of 
Korhonen et al.[p007) l and from this paper (see Fig. [5} are shown with 
green dots/arcs, green dots in the middle of the circles indicate polar 
caps. Phases of prominences (see Fig. [2|, i.e., when the prominences 
cross the stellar disk are shown with black arrows outside the circles. 
Phases of increased Ho emission seen in Fig. [2] are marked with red 
dots/arcs inside the circles. 


The longer (~7 years) cycle dominates the Her variations and 
seems to be shifted by about 1-2 years with respect to the bright¬ 
ness in V. With the growing spottedness between 1999-2002 
the Her excess decrease and spots dominate. From 2002 to 2004 
the star brightens which means less spots, by 2004 the Her is 
stronger, and after that more spots appear and the star becomes 
gradually fainter again. This could mean, that we see the same 
quasi-periodic behaviour in the photometry and the Her spec- 
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Fig. 8. Comparison between the area of Ha residuals from Gaussian fits 
and light curves. 
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Fig. 9. Magnitude-colour index curves of FK Com using photometric 
data of |Hackman et al.| < |2013| ) covering the years 1997-2010. Dashed 
lines show the best fit to the data. 





tra, i.e., the photosphere and the chromosphere, but the chro¬ 
mospheric changes are delayed with respect to the photosphere. 

|Radick et al.| ( |1998| ) investigated the long-term relations be¬ 
tween the photospheric and chromospheric variations of 35 stars, 
including the Sun, using the Mount Wilson HK data (see e.g. 
|Wilson|[l978| > with contemporaneous photometry from Lowell 
Observatory during the last decade until 1995. The results show 
that the younger stars of the sample, i.e., below 1-2 Gyr age are 
fainter when their chromospheric emission increase while the 
older stars, including the Sun, show parallel behaviour show¬ 
ing brightening with higher chromospheric emission. Note, that 
these relations are not strict parallels or minor images but show 
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Fig. 10. Light curve and magnitude-colour index curves of FK Com 
using photometric data of |Hackman et al7| < |2013} t covering the years 
1997-2010. The magnitude-colour index plots show four section of the 
light curve having lower amplitudes plotted over all data, marked with 
small dots. Dashed lines show the best fit to the given sections. 



Fig. 11. Top: AV photometry of FK Comae, bottom: Her line areas of 
the residual spectra. Continuous lines in the upper two plots show a 
two-period fit to the data (see text). 


Fig. 14. The upper three panels show the long-term behaviour of the 
b-y, B- V and V-Iq colour indices which reflect the global temperature 
changes of the star. In the bottom panel the starspot temperatures using 
photometry and Doppler Imaging are plotted. The data points for the 
Doppler-map temperatures are from this paper, Korho nen et al.| ( |T999| >, 
Korh onen et al.| < f2000) >, and |Korhonen et aTT| 4 2007 [ i . Horizontal line at 
5000 K shows the unspotted temperature. 


phase shifts and sometimes different patterns (see Radick et al. 
[19981 Fig- 3.). 

It is also possible that cool spots are caused not just by mag¬ 


netic activity but also vortices, as proposed by Kapyla et al. 
( |2011| l. In this case there would be no connection between the 
light curve and the Her emission. Another possibility could be 
the structure of the spots: i.e. spot groups with mixed polarities 
vs. large unipolar spots, as these two kind of structures could 
cause very different Her emission. In principle this could be ver¬ 
ified by studying the magnetic field structure. 

According to the Doppler images, bright regions are present 
on the surface. Such hot regions were found on FKCom in 
Doppler images of |Korhonen et al.) ( |2007|) , e.g. in 1993 June or 
2003 June, or in some of the ones presented in this paper (see 
Fig. [5}. These bright spots tend to get more pronounced when 
the data quality is poor, and especially large phase gaps in the 
data result in hot spots in Doppler images, suggesting they are 
artefacts of the inversion method. This possibility is also sup¬ 
ported by the fact that the bright regions mostly occur near the 
cool spots - this pairing of cool and bright features is a way for 
the inversion to mimic spectral features visible at certain phases 
but absent in others. This can also occur if there are significant 
changes in the spot configuration occurring during the observa¬ 
tion run (see Hackman et al .|2013[ ). Therefore, their reality can¬ 
not be conclusively verified. 

In Fig. 14 we plotted spot temperatures from both Doppler 
images and spot modelling. Photometric spot temperatures were 
derived both from narrow- and wide-band data ( b-y and 
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B-V,V-Ic respectively), using SpotModel (Ribarik et al. 2003[ i. 
The spot temperatures derived from the different colour indices 
generally match each other well. On average though, the spot 
temperatures from wide-band photometry are a little bit higher 
than from the narrow band data, which is probably due to more 
contribution of faculae to the wide-bands. The temperatures of 
the coolest areas from Doppler imaging agree fairly well with 
spot temperatures from photometry. The hottest regions have 
about 500 K higher temperatures than the photosphere (72ir = 
5000 K assumed). By 2010 the temperature of the coolest re¬ 
gions from Doppler imaging became higher supporting the mea¬ 
sured increase (blueing) of the V — Ic colour index. 

The contrast of the Doppler images seems to decrease in time 
in the good-quality maps. In Fig. 14 minimum and maximum 


going on in the activity of FK Com. This change might be con¬ 
nected with the results of Ko rhonen et al . (2009) and |Puzin et al.| 
( 2014| l, who studied longitudinal magnetic field on FKCom us¬ 
ing low resolution spectropolarimetry. According to Korhonen] 


et al. 2007). A linear fit to the coolest Doppler imaging temper¬ 


atures in these maps shows that the spot temperature increases 
on average approximately 23 K per year during the time pe¬ 
riod 1994—2010. On the other hand, the hottest temperatures in 
the maps are typically around 5600 K (67% of the time 5500- 
5700 K). The fit to the hottest temperatures exhibits a slight de¬ 
crease over the years, but only at a level of approximately 4 K 
per year, and it cannot be considered significant. One has to note, 
though, that the unspotted surface of FK Com has a temperature 
of approximately 5000 K, and therefore the hottest temperatures 
are significantly above this. 

The 2009/2010 observations are especially puzzling. The 
whole brightness level is low, the light curve has low amplitude 
modulation. As photometric observations show, a sudden phase 
jump or flip-flop occurred in the rotational modulation amount¬ 
ing about 130° in longitude, sometime in 2009 March. The phase 
of the main minimum jumped from phase 0.6 to 0.3 as seen from 
Fig 6. The Her spectral region, observed overlapping the photo¬ 
metric variations from late 2009/early 2010 shows very low level 
of chromospheric activity (see Fig. [8]). In Fig. [2] we can see that 
the level of Ho- emission is indeed steady, it is originating from 
outside the stellar disk, and there is a wave-like structure indi¬ 
cating a prominence, that reaches the middle of the stellar disk 
at phase 0.2-0.4, coinciding with the light curve minimum. 

In the end of the 2009/2010 season by 2010 April, the am¬ 
plitude of the rotational modulation became even lower and the 
main minimum shifted to phase 0.1. The spot temperatures from 
V - Ic was steadily decreasing from 2005 to 2010. In 2009 Au¬ 
gust the lower spot temperature from Doppler imaging (DI), ob¬ 
served later than the phase jump (or flip-flop?) happened in 2009 
March, is significantly warmer than that from the V - Ic colour 
index. Similarly, in the second half of 1997 a flip-flop was ob¬ 
served by Korhonen et al. ( |2001[ ), and after that, in early 1998 the 
lower spot temperature from DI was also warmer than those from 
photometry (at that time we have temperatures from 3 colour in¬ 
dices). These two observations are suggestive for finding a cor¬ 
relation between flip-flop events and chromospheric activity on 
a long-term base. 

To investigate further the relation of the bright and cool re¬ 
gions on FKCom, we compared magnitude-colour diagrams 
from different epochs. Four of these diagrams are shown in Fig. 
[TO] The V vs. V-Ic plots have very similar trends, but in the last 
two epochs this trend seems to be somewhat weaker (especially 
in the cases of lower amplitudes, plotted with a darker shade 
and with a separate fit). This can suggest, that there is a change 


et al. 


(2009)1, the average longitudinal magnetic field strength of 


FK Com was in the order of 200 G in 2008 April, but|Puzin et al. 


map temperatures in a uniform sample of Doppler images are 
given. All the maps used for this plot are from SOFIN data and 
they have maximum phase gap of < 0.25, i.e., cases where strong 
artefacts are not expected (temperatures are collected from this 
paper, Korhonen et al. 1999 Korhonen et al. 12000] and|Korhonen| 


( |2014) l found that in 2012 May there was no detectable magnetic 
field on the star. The measurements from 2011 April-May show 
significant detections of the longitudinal magnetic field, but with 
much reduced field strength compared to 2008 April (Korhonen 
et al. in preparation). The observations of |Puzin et al. (2014 1 are 
of lower accuracy than the ones used by Korhonen et al. (2009), 
and therefore, as the authors also hypothesise in their paper, the 
non-detection could be due to the field being too weak for them 
to detect. All this supports that the activity of FKCom has signif¬ 
icantly decreased during the recent years: decreasing spot con¬ 
trast of the Doppler maps, weak Ha levels in the last epochs, and 
decreasing longitudinal magnetic field strength. 

Describing the complex phenomena observed on FK Com is 
not an easy task. According to the most likely explanation, the 
object is a coalesced contact binary system, possibly with a sur¬ 
rounding disk ( ]Webbink||1976| , and often one or more promi¬ 
nence arcs. The stellar surface is covered both by bright plage 
regions, and dark starspots, as seen in the Doppler images. The 
spots in some cases can form a belt-like structure, resulting in an 
unvarying light curve that is fainter than the brightest state of the 
star (e.g. in 1993.31, see Korhonen et al. |2002 ) There is a con¬ 
nection between the changes in the chromosphere and the pho¬ 
tosphere, but it is far from obvious. According to the later light 
curve amplitudes, polarimetric measurements spot contrast from 
Doppler images, the general activity level seems to decrease, but 
the average light curve level still indicates a spotted state. 


5. Conclusions 

Based on the high and low resolution spectra presented in this 
work we draw the following conclusions: 

- The Ha region shows a double-peaked emission, possibly 
indicating a circumstellar disc, and often prominence-like 
structures that reach to more than a stellar radius and are sta¬ 
ble for weeks. 

- The low-resolution Ha spectra show small-scale fast vari¬ 
ability during one observing night on the time scale of hours, 
and the shape of the line profile changes between the nights 
more significantly. The first observing night, showing higher 
Ha emission level and a more asymmetric shape, coincided 
with the most spotted phase of the light curve and of the 
Doppler map; 

- Comparing the phased AT light curves and the excess Ha 
emission, we often find an anti-correlated behaviour of the 
two: the Ha emission peaks at the faintest state of the star: 
at the most spotted phase on the light curve and the Doppler 
maps; 

- This trend is not obvious on longer term: between 2000 and 
2006 the AV light curve and the Ha emission seems to be 
correlated, but the observations before and after this time 
seem to behave differently; 

- The long-term correlation between the light curve and the 
Ha emission could be explained by the plage regions, but 
the magnitude-colour index in Korhonen et al.| ( |2001| ) does 
not confirm this scenario: the star seems to get redder as it 
gets fainter, indicating a spot-dominated activity; 
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- The connection between the Her emission level and the light 
curve could be explained by the same quasi-periodic changes 
in the chromosphere and the photosphere, with the chromo¬ 
spheric changes being 1-2 years Tate’ compared to the pho¬ 
tosphere; 

- In the 2009/2010 observations the average level of the AV 
light curve is very low compared to low level of the Her emis¬ 
sion and the low contrast of the Doppler maps; 

- On the whole, the activity of FK Com is unusually weak in 
2009-2010. 
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Appendix A: Comparison between the different 
activity features 

In the following, we give a detailed description of the Her activity 
in those cases where indications of a prominence was seen, and 
their possible connection with photospheric features: 


1998.27/1998.30 The prominence crossed the stellar disk at » 
0.8 - 0.9 phase. The 1998.27 set shows steady Ha emission, 
but there is a flare in the 1998.30 set at phase 0.9. The light 
curve shows a minimum around 0.75 phase. The Doppler im¬ 
age in Korhonen et al. ( 2007| show high-latitude spots at 
phases 0.34-0.52 and 0.78-1.00. The flare, the dark spots 
and the prominence coincide in phase. 

1999.41/1999.58 In the 1999.41 observations we see a promi¬ 
nence at phase « 0.2. The light curve shows a minimum 
at phase 0.9 and a maximum at 0.4. The Doppler image 
from 1999 May shows a polar spot, and high-latitude spots 
at phases 0.6 and 1.0 and a southern spot at 0.97, but the one 
from 1999 July shows a much weaker spottedness around 
the same phases, and no polar spot. In this case the promi¬ 
nence is close to the brightest phase of the star. In both May 
and July there are brighter regions in the maps around 0.25 
phase: in May we see two bright areas on the northern hemi¬ 
sphere, while in July it was found on the southern part. 
1999.58 is the only case when we see a strong indication of 
two coinciding prominence loops in the Her dynamic spectra, 
one around 0.5-0.6 phase and one at phase 1.0. 

2000.33 We see a prominence crossing the stellar disk at 0.7 
in phase. The light curve at this epoch shows a minimum 
at phase 0.7 and a maximum at 0.3. Doppler map from 2000 
March shows high-latitude spots around phases 0.5-0.64, the 
one from 2000 April has dark features around 0.23, 0.52- 
0.66, 0.59-0.70,0.78-0.84, 0.97 phases in different latitudes. 
This latter map was marked in the paper as less reliable. The 
spotted regions between phases 0.5 and 0.8 are present in 
both maps, one of these might be connected with the promi¬ 
nence. 
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2001.35 The light curve in this epoch shows a very symmetric 
sinusoidal shape with a minimum at 0.46 and a maximum 
at 0.97 in phase. In this case the wave-like Her structure at 
phase 0.5 in Fig. [2] is less obvious, probably the promi¬ 
nence is hidden by the stellar disk around phase 1.0. This 
suggests a low-latitude feature. The Doppler map from 2001 
May shows the strongest features - both dark and bright re¬ 
gions - around 0.25 phase, although this map was marked as 
less reliable. The following map from 2001 June shows two 
high-latitude spots at phases 0.20-0.38 and 0.54-0.74, an 
equatorial spot at 0.49, and a brighter feature around phase 
0.5. 

2004.10 The light curve in this epoch is a sinusoidal one with 
a minimum at 0.95 and a maximum at 0.43. The closest 
Doppler image is from 2004 July, shown in Fig.[5](in July the 
shape of the light curve is similar to that in February, but the 
amplitude is higher and the maximum is around phase 0.55). 
In this map the most pronounced features are two dark spots: 
a mid-latitude one around phase 0.1 and a high-latitude one 
between phases 0.2 and 0.3. The dynamic spectrum shows 
that the prominence is crossing the stellar disk at ~ 0.5 phase. 
The wave-like shape in the radial velocity is somewhat dis¬ 
torted, which could be a result of a mass ejection. We can see 
an increased level of Her emission at phase 0.7, that could 
be originating from the prominence. According to both the 
Doppler maps and the light curves in this case there is possi¬ 
bly no connection between the dark photospheric spots and 
the prominence. 

2004.58 The prominence crosses the stellar disk around phase 
0.05, that is close in phase both to the light curve minimum 
and the large dark spots seen in the Doppler maps. 

2007.56 The light curve shows a minimum at phase 0.8 and a 
maximum around phase 0.27. The Doppler map in Fig. [5] is 
very interesting, it is dominated by both bright and dark spots 
between phases 0.55-1.00, with dark spots from the southern 
hemisphere reaching to higher latitudes, and a bright elon¬ 
gated feature at ~ 45° latitude. The dynamic spectrum shows 
a broad wave-like feature, indicating a prominence crossing 
the stellar disk around phase 0.6, close both to the light curve 
minimum and the spots found in the Doppler maps. 

2009.67 The Doppler maps from 2009.67 and 2010.00 are quite 
similar: both have a spot group around phase 0.25 at rela¬ 
tively low latitudes. The temperature contrast is low, about 
500 K cooler than the unspotted surface in 2009 August and 
600 K as in 2010 January. There is no sign of a polar spot. 
There is no photometric observation during this 2009.67 pe¬ 
riod, but the light curve before these observations is sinu¬ 
soidal one with a minimum at phase 0.32 and a maximum 
at 0.78. The following AV data shows similar behaviour as 
the previous one, with an additional secondary minimum 
around phase 0.65. The dynamic spectrum is very interest¬ 
ing: it shows an obvious wave-like structure, but the general 
level of the Her emission is very low, barely higher than the 
lowest observed one. The source of the emission crosses the 
stellar disk at phase 0.3. 
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Fig. 13. AV light curve (top) and its time-frequency analysis using the 
STFT (Short-time Fourier transform) method (bottom) showing multi¬ 
ple activity cycles of changing length on FK Comae. 


Fig. 12. Fourier spectrum of the FK Com V light curve in those 
year, where Her observations exist, a) shows the whole spectrum, b) 
is zoomed in to the long-period region, c) shows the zoomed-in spec¬ 
trum after pre-whitening with the 6.8 year-long period. The bottom plot 
d) shows the spectral window. 






































